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We reexamine the effects of long-range Coulomb interaction on the onset of superconductivity. 
We use the model of N complex scalar fields with the Coulomb interaction studied first by Fisher 
and Grinstein (FG). We find that near d=3 space dimension, the system undergoes second order 
phase transition if TV > 55.39, but undergoes possible fluctuation driven first order transition if 
N < 55.39. We give the detailed derivation of the field theory Renormalization Group (RG) of this 
model to one loop. Our RG results disagree with those of FG near d — 3. Possible scenario at d = 2 
is proposed. 



I. INTRODUCTION 



Superconductor to insulator (SI) transitions have been observed in disordered thin films by systematically varying 
the thickness of the filmsu or by tuning the magnetic fieldsB. A closely related SI transition has also been studied 
in artificially constructed 2d Josephson-junction arrays by systematically varying the ratio of charging energy and 
Josephson coupling energyu or by tuning the magnetic fieldsu. Universal conductivities have been measured right at 
the SI transitions in theseusystems. More recently, SI transition was demonstrated in 3d granular Al-Ge samplesa. 
Many theoretical effortsH~tL3 have been made to investigate the SI transitions in these experimental systems. It is 
generally accepted that the transition is correctly described by a model of interacting charge-2e bosons moving in a 2d 
random potential. A related, but much simpler model is Boson-Hubbard model which consists of bosons hopping on a 
periodic lattice. When the boson density is commensurate with the lattice, this model displays a quantum transition 
from Mott insulator to a superconductor. Neglecting disorder as-a-Stst approximation, a lot of authors studied the SI 
transition and its universal conductivity in this simplified model0O~cij. Near the quantum critical point, the effective 
low energy theory of this model is simply a </> 4 theory: 
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where <f> m are m = 1 • • • N species of charge e* = 2e complex bosons. 

In order to confront with experiments, several important effects such as disorder, dissipation and long-range Coulomb 
interaction should be taken into account. From general scaling arguments, Fisher, Grinstein and Girvin argued that 
the long range. Coulomb interaction between the charged bosons will render the dynamic exponent z to be 1 
any dimension!!! This is indeed what was observed on the magnetic field tuned SI transition in thin 2D film: 
The measurements on magnetic field-induced SI transitions in Josephson- Junction arrays are also consistent with 
z = 1, v > 1. Therefore, it is very important to investigate the effects of Coulomb interaction on SI transition in any 
detail. In Ref.0, Fisher and Grinstein (FG) studied the model of charged bosons with Coulomb interaction hopping 
on a lattice. By taking continuum limit, they found that the low energy effective theory is described by 
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Where ^4o is the time component of the U(l) gauge fields and is responsible for the long-range interaction between 
the charged bosons. The X/r Coulomb interaction corresponds to a = d X. ,. . 

The model Eq.|| should not be confused with the model studied in RefsJli'taO. The crucial difference is that only 
the time component Aq of gauge fields is i nvolv ed in Eq.|^, therefore, the Lorentz invariance is broken. However, in 
the scalar electrodynamics studied in Refs.c3u3't2l, both time and transverse components are involved, the Lorentz 
invariance is respected. 

In principle, the coupling of the order parameter to the spatial components of the U(X) gauge fields should also be 

A 2 

included in the above equation. However, at 2D films, the penetration depth \2d — with A the bulk penetration 

n o o o 

depth and d the film thickness. In the experimental systems^, A ~ 100 A, d ~ 5 A, so \2d ~ 2000 A, the spatial 
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coupling will not manifest itself until experimentally unobservably close to tha SI transitionBilll. Therefore, it is safe 
to neglect the spatial coupling in Eq.|2[ In 3D samples investigated in Reffl, we assume that they are in extreme 
type-II limit. In this limit, the spatial coupling can also be neglected. 

FG did Renormalization Group (RG) analysis of the model by performing a double expansion in e = 3 — d and 
e a = 2 — a. The 1/r Coulomb interaction corresponds to e = e^. They concluded that depending on parameters, 
at two space dimension, the transition can be either first order or second order with Coulomb coupling marginally 
irrelevant at the 3D XY critical point. 

In this paper, we provide a detailed RG analysis of the Model Eq.|| using field theory method. We do not perform 
a double expansion in e and e CT . Instead, for simplicity, we fix a to be 2 ( namely e CT = ) and perform a e = 3 — d 
expansion, a = 2 corresponds to 1/r Coulomb interaction at d = 3 and logarithmic instead of 1/r interaction at 
d = 2. We find that near 3 space dimension, the system undergoes second order phase transition if N > 55.39, but 
undergoes possible fluctuation-driven first order transition if N < 55.39. We calculate the dynamic exponent z, the 
boson propagator exponent 77, the gauge field ( Coulomb interaction ) exponent tja and the correlation exponent v. 
We give the detailed derivation of field theory RG of this model, because the RG method we used is interesting in its 
own right. We expect the structure brought out by our RG procedpses have some general impacts on how to perform 
correct RG on other zero temperature quantum critical phenomenal!^. The same method has been applied successfully 
to study the quantum transition from Fractional Quantum Hall state to insulating state in a periodic potentials. 

If putting £p to zero in Eq. (7) of Ref.Q, we find our RG results disagree with those of FG. Although no details are 
given in Ref.Q, we suspect that (1) the anisotropy between space and time (namely the dynamic exponent z ) was not 
treated correctly by FG (2) The Ward identity ( Eq^ in Sec. II ) is violated in the momentum shell method employed 
by FG. 

In this paper, we do not even intend to perform the double expansion in e = 3 — d and e CT = 2 — a. Because we 
believe that when gauge field fluctuations are involved, extrapolating to the physical case e = e a — 1 can only lead to 
misleading results. 

The rest of the paper is organized as follows. In Sec. II, we introduce the model and set up the formulation to 
perform RG. In Sec. Ill, we give the detailed RG analysis and bring out its elegant structure. We find that our results 
disagree with those of FG and point out the possible mistakes made by FG. In Sec. IV, we give more discussions and 
point out the possible scenario at d = 2 and some future directions. Finally, in the appendix, we perform RG on the 
well-studied Scalar Electrodynamics in Feymann gauge as a non-trivial check on the correctness of our calculations 
in Sec.III. 

II. THE MODEL AND FORMULATION 

In order to perform RG analysis, we rewrite Eq.|^ in the following formEl 
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Where Xi (r) are spatial (temporal) co-ordinates with do = d T , di = d Xi . we are working in d = 3 + e spatial 
dimensions and [i is a rcuorpalization scale. The parameter a is introduced to allow for anisotropic renormalization 
between space and timeElrc3. 

It is easy to see Eq.^| is invariant under the space-independent gauge-transformation 

-> cjje^ , A Q A Q - ^doA(r) (4) 
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In the perturbation theory, the three kinds of vertices in Fig.l are needed 
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FIG. 1. The Feymann rules for three kinds of vertices in Eq.H (a) the type 1 vertex (b) the type 2 (seagull) vertex (c) the 
quartic coupling of bosons 

The loop expansion requires counterterms to account for ultraviolet divergences in momentum integrals; we write 
the counter terms as 
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The Ward identity following from the gauge invariance dictates 

Z\ = Z a 

Using the identity, we relate the bare fields and couplings in S to the renormalized quantities by 

4>mB = Zl /2 (j) m 
,1/2 



(5) 



(6) 



A b = Z 3 A 
as = {Z a /Z 2 ) 1/2 a 
ef i e/2 (Z a /Z 2 ) 



1/4^-1/2 



e 7 ^—3/2^-1/2 

ub = u[i Z^Z 2 Z a 1 
The dynamic critical exponent, z is related to the renormalization of a byfi 

i d 1 1 1 d 1 22 
z = 1 — fJ,— m a = 1 i^—r- l n — 
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We define the exponent r\ by the equal-time correlation function 

< 0(x,r)0 t (O,r) >~ x -(d+z-2+v) 
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It is easy to see this correlation function is invariant under the gauge transformation EqJJ, therefore, r/ is indeed a 
gauge invariant quantity and given by 



V = A t_ T" m -^2 

The anomalous dimension t\a of the gauge field is also a gauge invariant quantity: 



(10) 
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Near the critical point, the gauge field propagator is 

< A (-k,-uj)A (k,u)) >r 
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Physically, it means that at the critical point, at d spatial dimension, the long-range interaction takes the form 

-(d-2+r, A )_ 

Finally, the critical exponent v is related to the anomalous dimension of the composite operator <jr m <p m by 

^-2 = /^ln^ (13) 
The renormalization constant can be calculated by inserting the operator into the boson self-energy diagrams. 

III. THE PERTURB ATIVE RG CALCULATION AT ONE-LOOP 
A. The boson self-energy 

First, we consider the boson self-energy diagrams in Fig. 2 
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FIG. 2. The self-energy diagram of bosons 



Using the Feymann rules listed in Fig.J^, we can bring out the divergent structure of Fig. 2a by conventional Dimen- 
sional Regulization (DR) method 
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Where • • • means all the finite terms. 

In the above equation, we scaled all the frequencies by the anisotropy parameter a, so a does not appear explicitly 
in the above and in the following. 



Similarly, we can evaluate Fig. 2b and 2c 
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In the last equation, the convention of DR is used. 
Adding all the contributions from FigJ|, we get 



(2 fl ) + (26) + (2 C ) = i - 5 -(a; 2 - y ) + . 
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It is easy to see that we encounter the UV frequency divergences in both Fig. 2a and Fig. 2b. However, they carry 
opposite signs, therefore, cancel each other. Similar cancellations also appear in the Dirac fermion model studied in 
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Rcf.0 and will appear in the following calculations of other Feymann diagrams. We expect this cancellation of UV 
divergences in frequency is a general feature of non-relativistic quantum field theories describing zero-temperature 
quantum phase transitions. 

From Eq.^, we can identify the two constants Z a , Z 2 in Eq.|^ in the Minimal Subtraction scheme 

e 2 e 2 

z « = 1 + t^> Z2 = 1 -7Tir ( 17 ) 
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Note Z a ^ Z 2 . This is due to the lack of Lorentz invariance of our model. 



B. The type 1 vertex 

We turn to the evaluation of the renormalization of type 1 vertex in Fig. 3 
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FIG. 3. The one-loop diagrams of type 1 vertex 



Applying the Feymann rules to Fig. 3a, we have 
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It is easy to see Fig. 3b vanishes identically and Fig. 3c and Fig. 3d are 
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(3c) + (3d) 
Overall, we find for the type 1 vertex 
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(3a) + (36) + (3c) + (3d) 
The constant Z\ in Eq.|| can be identified 

Zx = l- 



4tt 2 . 
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From Eqs.[l7],|2l], it is evident that the Ward identity EqJ^ indeed holds. The Ward identity can also be shown by 
evaluating the renormalization of the seagull vertex. 



C. The polarization of gauge field 

Next, we calculate the vacuum polarization graph of the U(l) gauge field in Fig. 4 

k+q 





(a) (b) 
FIG. 4. The vacuum polarization of gauge fields 

The fermion bubble has Lorentz invariance, therefore, the momenta in Fig. 4 are D = d + 1 = 4 — e dimensional 
momenta 
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where N comes from the summation over the boson suffix m. 

In Eqj^, only the time component of the gauge field is involved. Putting /i = v = in the above Eq. we find 
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where k is d dimensional space momentum. 
The constant Z% in Eq.|| can be identified as 
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D. The renormalization of the quartic term 

In this subsection, we compute the renormalization of the quartic term. 
First, let us look at the contributions from the quartic coupling in Fig. 5. 
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FIG. 5. The renormalization from the quartic coupling 

As in FigJ^, the momenta in the above figure are also D dimensional momenta. 

, . u 2 u 2e f d D q 1 
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Similarly, we can calculate the contributions from Fig. 5b and 5c@ 
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(56) = (5c) = 
Adding all the contributions from FigJ|, we get 
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Second, let us look at the contributions from the seagull term in Fig. 6 
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FIG. 6. The renormalization from the seagull term 
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As in Eq.|TJ and we encounter the UV divergences in the frequency integrals. We leave them alone at this 
moment and go ahead to calculate the diagrams in Fig. 7 





(a) (b) 
FIG. 7. The renormalization from the combination of the type 1 term + the quartic term 



The corresponding expressions are 
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It is important to note that Fig. 7a and Fig. 7b have opposite signs. Actually, there are two diagrams in class (7a) 
and four diagrams in class (7b) corresponding to different ways to put the photon line, so the overall contributions 
are 
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Let us look at the contributions from the type I vertex in Fig 
1 3 1 
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FIG. 8. The renormalization from the type 1 term 
The corresponding expressions are 

(a ) — ( f ^9 f dv ( V + ^lX^ + ^3X^1 + 2^2 ~ V)(UJ! + UJ 2 + tOj ~ v) 
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Actually, there are two diagrams in class (8a) and two diagrams in class (8b) corresponding to different ways to 
put the photon line, so the overall contributions from Fig. 8 are 



2(8a) + 2(86) = 2{e^ 2 ) A 



d d q 



-F7T(- 



1 



1 f dv 3e 4 

(2^) d ( ( j'_fc 1 )2 V (^_fc 3 )2 ' (q - fc 4 )2 7 27 _ 2^ + ' ' 

Finally, we need to compute the diagrams in Fig. 9 
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FIG. 9. The renormalization from the type 1 term + the seagull term 



The corresponding expressions are 
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Exchanging the leg 3 with leg 4 of Fig. 9a and 9b, we get another two diagrams (9c) and (9d) which were not shown 
explicitly 
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The over all contributions are 
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Adding all the contributions from Eqs. 
cancel and lead to the constant in Eq 




we find the UV divergences in the frequency integrals indeed 
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E. The calculation of (3 function 

In the previous subsections, we did the explicit calculations of the renormalization constants by considering a direct 
perturbative expansion in the Coulomb fine structure constant w = e 2 and the quartic coupling u. At one loop order, 
the values of the renormalization constants are summarized as 

in 

z a = 1 
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Substituting the above equation into Eq.[?], we find the (3 functions for w and 
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Where A = VA 2 - 52A - 188. 

It is easy to see that A is imaginary for N < N c = 55.39. This A c is much smaller than the corresponding critical 
number n c = 365 in the Scalar Electrodynamics studied in Rcf.EJ. For N < N c , u is complex for both fixed points and 
physically inaccessible. The only accessible fixed points are the Gaussian and Heisenberg fixed points, both of which 
are unstable to turning on Coulomb interaction. There is a runaway flow to negative value of u which was interpreted 
as a fluctuation-driven first order transition in RefJla. For A > A c , these fixed points have real values for u* and are 
physically accessible. 

In order to compare our (3 functions with those of FG, we put A = 1 in Eq.|38| 

5 

(3( w ) = -ew + —w 2 

„, . 5u 2 wu w 2 , . 

^=- eU+ T^-4^ + 2^ (40) 

Setting w — ^-Wj g ,u — 32iT 2 Uf g , the above equations become 

1 



P( w fg) = -ew/s + 3 W / 9 



2 1 

P( u fg) = ~ eu fg + 10u /s ~ l w f9 u fa + ( 41 ) 

Putting e CT = in Eq.(7) ofRcfi, we find the R. G. flow equations got by FG are different from ours. The j3 function 
should be gauge- independent at least at one- loop order ( see also the discussions in the appendix on this point). No 
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details being given by FG, we can only guess the possible mistakes made by FG: (1) the anisotropy between space and 
time (namely the dynamic exponent z ) was not treated correctly by FG (2) The Ward identity in the time component 
Eq.|| is violated in the momentum shell method employed by FG. 

The disagreement at e a = puts doubts of the conclusions reached by FG at Coulomb interaction case e = e a . 

Let us look at the solution right at the critical dimension d = 3, Putting e = in Eq.38 and taking the ratio of the 
two equations, we find the trajectories of the RG flow satisfy the homogeneous differential equation: 



du 3,u. 2 6m 12 
dw ~ 2w N + 4w + N + 



(42) 



If N > N c , the solution of this equation is 
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Where C is an arbitrary constant. All the trajectories flow to the origin with two marginally irrelevant couplings w 
and u. At finite temperature, they will lead to logarithmic corrections to naive scaling functions 
If N < N c , the solution of this equation is 



- . 2(A^ + 4) r 3(A^ + 4).M 2 1 

C + In«, = arctan[A_^ ( _ _ (__ + _))] 

All the trajectories flow to the negative values of u, indicating a fluctuation driven first order transition. 



(44) 



F. The calculation of critical exponents 

In this subsection, we calculate the critical exponents of the second order transition when N > N c 
From EqsMpAp9L we find the dynamic exponent z 



z = l- 



4c 



From Eqs.|T^,^,^9|, we find the exponent t] 



From Eqs]ll]|37y39L we find the exponent t]a 
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(47) 



The long-range interaction is modified from the form ~ r 1+6 to the form ~ r 1 + 4e /( JV + 4 ) p 
Z^di can be calculated from Fig. 10 





(k,C0) (q,v) (q,v) 



(a) (b) 

FIG. 10. The cross stands for the insertion of the composite operator at zero momentum 
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The corresponding expressions are 
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where we put the mass of the boson m into the boson propagators in order to avoid the infrared divergence. 
In contrast to the self energy diagrams Fig. 2, there is no UV frequency divergence in the above equations. 
From the above equation, we can identify 

u(N + l) e 2 , . 

Substituting the value of Z2 in EqjlT], we get 



_ u(N + l) w 

+ (50) 



The correlation length exponent v is given by 



t , N+ 1 r N - 1 6 , A 



v~ = 2 ; -e\ + ± 1 (51) 

2(7V + 4) l JV+l JV + 4 JV + 4 1 v ' 



IV. DISCUSSIONS AND CONCLUSIONS 



In this paper, we provide a rather detailed derivation of RG analysis of Eq.|| We find the delicate cancellation of the 
UV frequency integral divergences, therefore, we only need to regulize the UV divergences in momenta integral. We 
believe that this is due to the space-independent gauge invariance Eq. | ( or the Ward ideatity in the time component 
Eq.|| ). Similar cancellations also appear in the Dirac fermion mode studied in RefEJca. We expect this is a very 
general feature of zero temperature quantum critical phenomenally. This feature in non-relativistic quantum field 
theory reminds us of the theorem in relativistic quantum field theory that infra-red divergences must be cancelled out 
in any physical processes. 

Here, we would like to mak&-a_brief comparison between the boson model studied in this paper and the Dirac 
fermion model studied in Refcj'ca. In Dirac fermion model, both Coulomb and Chern-Simon couplings become 
marginal at d = 2, the competition between this two couplings produce a line of fixed points with z = 1 and non- 
vanishing renormalized Coulomb coupling. The fermion quartic term is irrelevant at the non-interacting fixed point 
and remains irrelevant on this z = 1 line of fixed points. In the boson model studied in this paper, both the Coulomb 
and the boson quartic couplings become marginal at d = 3. Fixing a to be 2 ( namely e CT = ), we performed a 
e = 3 — d expansion, a = 2 corresponds to 1/r Coulomb interaction at d = 3, l/r d ~ 2 interaction at space dimension d 
and logarithmic instead of 1/r interaction at d — 2. The Coulomb coupling was found to drive the quartic coupling 
to negative value, indicating a possibls-&rst order transition. 

Halperin, Lubensky and Ma (HLM)Hj investigated the model describing the Superconducting to Normal (SN) and 
Nematic to Smectic-A transition in d = 3 bulk system. By performing e = 4 — d expansion, they found a new stable 
fixed point with non- vanishing charge if the number of the order parameter components N > M c — 365.9/2, but 
runaway RG if TV < M c (see also the appendix). HLM interpreted this runaway RG flow as indicatingJiuctuation- 
driven first-order transition. Later, Dasgupta and Halperin(DH) stuided the model directly on a 3d latticed. By using 
Monte carlo and duality arguments, they found the transition should be a second order transition in the universality 
class of the inverted XY-model instead of a first-order transition. DH's results indicate that the e expansion may 
break down at d = 3. 

In this paper, we find the critical value _/V c in the pure Coulomb case is much smaller than the corresponding 
value M c in the 3d SN transition ( more precisely AD Scalar Electrodynamics ). This shows that the fluctuation is 
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considerably reduced in the pure Coulomb case. We believe that it is more likely the SI transition at d — 2 should be 
governed by a non-trivial fixed point with non-vanishing Coulomb coupling and dynamic exponent z = 1 

In order to understand the effects of 1/r interaction at d — 2 near the 3D XY fixed point, two possible expansions 
can be used. One route is to fix at d — 2 and perform large TV expansion. Unfortunately, as shown in this paper, the 
direct large N expansion can only probe the physics at N > AL— -therefore, is not very useful at physical case N = 1 in 
sharp contrast to the Dirac fermion case as discussed in Ref.Effia. Another route is the double expansion in e = 3 — d 
and £ a = 2 — a performed by FG. However, we believe that when gauge field fluctuations are involved, extrapolating 
to the physical case e = e a = 1 can only lead to misleading results- A. method directly works at D = 2 + 1 is needed. 

In order tojnake serious attempts to compare with experimentsoij, disorder has to be incorporated into EqJ^. As 
shown in Ref.El, due to Griffith effects, even weak disorder presumably produces a gapless Bose glass phase between the 
Mott insulator and superconducting phase. It is certainly very important to study the effects of Coulomb interaction 
on the transition from superconductor to this boson glass. It was suggested that with Coulomb interaction, the 
dynamic exponent z should be equal to 1 at all d and the renormalized Coulomb coupling is finite at the-transitionQ. 
This scenario was confirmed by extensive Monte-Carlo simulations of bosons in a 2d disordered mediumEZl. It will be 
welcome if we can establish this scenario analytically on this concrete model. This scenario was indeed, established in 
a clean 2d lattice model displaying quantum transitions from Quantum Hall state to insulating statecJ. 
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APPENDIX A: THE CALCULATION OF SCALAR ELECTRODYNAMICS IN FEYMANN GAUGE 

The Lagrangian studied for Scalar Electrodynamics in RefslllBIll is 

C = - ieA^ml 2 + +r\4> m \ 2 + J/i e |</> m | 4 



\(F,.uf + ^{d,A,f (All 



where Fu V = d^Av - d^A^. 

EqJAy describes transitions from superconductor to normal metal and from a smectic-A to a nematic in liquid 
crystals at d = 3. 

In all the previous work, the calculations were done in Landau gauge where A = 0. The advantage of The Landau 
gauge is that the divergent parts of Figs. 7, 8, 9 all vanish. In this appendix, we do all the calculations in Feymann 
gauge where A = 1. In this gauge, Figs. 7, 8, 9 do make contributions, but the geometrical factors in front of all the 
integrals should be the same as those in the Coulomb interaction case, therefore, by comparing our final results with 
those obtained in Landau gauge, we can make a non-trivial check on our RG results. Our final results show that (1) 
the (3 function is indeed gauge independent at least to one-loop (2) the calculations done in the Coulomb interaction 
case is correct. 

The loop expansion requires counterterms to account for ultraviolet divergences in momentum integrals; we write 
the counter terms as 

C C T = (Z 2 - l)|^ m | 2 + ^(Z 4 - l)fi e a\<j> m \ 4 

- i{Z x - l)e^ 2 A^d^l0m - 4>\ n dM + (Zi - l)e 2 ^Al4 n cb m 

+ i(Z 3 -l)(^) 2 (A2) 
Lo rcnt z invariance requires Z a = i?2, therefore, only one constant Z^ is needed. 



Eq.Al is invariant under the gauge-transformation 

0eiA(*) i A^A^- ^A(x) (A3) 

The Ward identity following from the gauge invariance dictates 

Zi = Z 2 (A4) 
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Using the identity, we relate the bare fields and couplings in C to the renormalized quantities by 



PmB 


7 l/2, 
— ^2 Vm 


Af,B 


- 7 1/2 A 


e-B 


= e^Z- 1 


UB 




X b 


= \z 3 



(A5) 



At one loop order, the values of the renormalization constants are 

w 



Zi = 

z 3 = 
z 4 = 
— 



1 + 

1 - 

1 + 
1 + 



47r 2 e 
Nw 



247r 2 e 
(N + 4)u 

167r 2 e 
(N+l)u 

16tt 2 c 



4ir 2 e 
3w 
8^~e 



3w 2 
2n 2 ue 



(A6) 



Substituting the above equation into Eq.A5, we find 

N 



(3(u) 



-ew - 



tW 



24tt 2 
(N + 4)u 2 

167T 2 



3wu 



3w 2 
2ir 2 



(A7) 



Eqs.A6,A7 should be compared with the corresponding Eqs.^,^8| in Coulomb interaction case. The differences 
should be noted. 



After proper scalings of w and u, we find EqA7 is the same as that derived in Landau gauge where A = 0. 
The critical exponents are 



6 

Vf = ( 

' N 

i]A = e 



y- 1 =2 



2(N + A) ly ' N N 



A„] 



(A8) 



Where A Q = VN 2 ~ 180A - 540. 

It is easy to see that r\ differs from the value computed in Landau gauge ( rj^ = — -^e ). It is not expected to be, 
because the correlation function of (f> is not gauge invariant anyway. What is the physical meaning of 77 is not clear. 
HowEKejj, the r\ in Eq.[l6] calculated in Coulomb interaction case is well-defined. The reader is encouraged to read 
Ref.ErEj for the similar discussions in Dirac fermion case. 

r]A, v are the same as those computed in Landau gauge, because they are gauge invariant quantities. r\A — e indicates 
that the Coulomb interaction remains 1 jt to order e at d -j . 3 i — e in contrast to the pure Coulomb interaction case. 
Actually, it was shown that the equation r\A — e is exacEQ'cQ. The critical number which divides the second order 
phase transition and the possible fluctuation-driven first order transition is M c = 365.9/2 which is of course gauge- 
independent. 
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